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EDITORS' PREFACE. 



The volume called 'Higher Mathematics, the first edition' 
of which was published in 1896, contained eleven chapters by 
eleven authors, each chapter being independent of the others, 
but all supposing the reader to have at least a mathematical 
training equivalent to that given in classical and engineering 
colleges. The publication of that volume is. now discontinued 
and the chapters are issued in Separate form. In these reissues 
it ' will generally be found that the monographs are enlarged 
by additional articles or appendices which either amplify the 
former presentation or record recent advances. This plan of 
publication has been arranged in order to meet the demand of 
teachers and the convenience of classes, but it is also thought 
that it may prove advantageous to readers in special lines of 
mathematical literature. 

It is the intention of the publishers and editors to add other 
monographs to the series from time to time, if the call for the 
same seems to warrant it. Among the topics which are under 
consideration are those of elliptic functions, the theory of num- 
bers, the group theory, the calculus of variations, and non- 
Euclidean geometry; possibly also monographs on branches of 
astronomy, mechanics, and matliematical physics may be included. 
It is the hope of the editors that this form of publication may 
tend to promote mathematical study and research over a wider 
field than that which the former volume has occupied. 

December, 1905. 
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AUTHOR'S PREFACE. 



This brief sketch of the Harmonic Functions and their use 
in Mathematical Physics was written as a chapter of Merriman 
and Woodward's Higher Mathematics. It was intended to give 
enough in the way of introduction and illustration to serve as 
a useful part of the equipment of the general mathematical 
student, and at the same time to point out to one specially inter- 
ested in the subject the way to carry on his study and reading 
toward a broad and detailed knowledge of its more difficult 
portions. 

Fourier's Series, Zonal Harmonics, and BessePs Functions of 
the order zero are treated at considerable length, with the inten- 
tion of enabling the reade'r to use them in actual work in physical 
problems, and to this end several valuable numerical tables 
are included in the text. 

Cambridge, Mass., December, 1905 
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HARMONIC FUNCTIOMS. 



Art. 1. History and Description. 

What is known as the Harmonic Analysis owed its origiiv 
and development to the study of concrete problems in various 
•branches of Mathematical Physics, which however all involved 
the treatment of partial dififerential equations of the same 
general form. 

The use of Trigonometric Series was first suggested by 
Daniel Bernouilli in 1753 in his researches on the musical 
vibrations of stretched elastic strings, although BesseFs Func- 
tions had been already (1732) employed by him and by Euler 
in dealing with the vibrations of a heavy string suspended from 
one end; and Zonal and Spherical Harmonics were introduced 
by Legendre and Laplace in 1782 in dealing with the attrac- 
tion of solids of revolution. 

The analysis was greatly advanced by Fourier in 18 12-1824 
in his remarkable work on the Conduction of Heat, and im- 
portant additions have been made by Lam^ (1839) ^"^ hy a. 
host of modern investigators. 

The differential equations treated in the problems which 
have just been enumerated are 
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^ HARMONIC FUNCTIONS. 

for the transverse vibrations of a musical string ; 

for small transverse vibrations of a uniform heavy string sus- 
pended from one end ; 

which is Laplace's equation ; and 

for the conduction of heat in a homogeneous solid. 

Of these Laplace's equation (3), and (4) of which (3) is a 
special case, are by far the most important, and we shall con- 
cern ourselves mainly with them in this chapter. As to their 
interest to engineers and physicists we quote from an article 
in The Electrician of Jan. 26, 1894', by Professor John Perry: 

" There is a well-known partial differential equation, which is 
the same in problems on heat-conduction, motion of fluids, the 
establishment of electrostatic or electromagnetic potential, certain 
motions of viscous fluid, certain kinds of strain and stress, currents 
in a conductor, vibrations of elastic solids, vibrations of flexible 
strings or elastic membranes, and innumerable other phenomena. 
The equation has always to be solved subject to certain boundary 
or limiting conditions, sometimes as to space and time, sometimes 
as to space alone, and we know that if we obtain any solution of a 
particular problem, then that is the true and only solution. Further- 
more, if a solution, say, of a heat-conduction problem is obtained 
by any person, that answer is at once applicable to analogous prob- 
lems in all the other departments of physics. Thus, if Lord Kel- 
vin draws for us the lines of flow in a simple vortex, he has drawn 
for us the lines of magnetic force about a circular current: if 
Lord Rayleigh calculates for us the resistance of the mouth of an 
organ-pipe, he has also determined the end effect of a bar of iron 
which is magnetized; when Mr. Oliver Heaviside shows his match- 
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HISTORY AND DESCRIPTION. 9 

less skill and familiarity with Bessel's functions in solving electro- 
magnetic problems, he is solving problems in heat-conductivity or 
the strains in prismatic shafts. How difficult it is to express exactly 
the distribution of strain in a twisted square shaft, for example, and 
yet how easy it is to understand thoroughly when one knows the 
perfect-fluid analogy! How easy, again, it is to imagine the electric 
current density everywhere in a conductor when transmitting alter- 
nating currents when we know Mr. Heaviside's viscous-fluid analogy, 
or even the heat-conduction analogy! 

" Much has been written about the correlation of the physical 
sciences; but when we observe how a young man who has worked 
almost altogether at heat problems suddenly shows himself ac- 
quainted with the most difficult investigations in other departments 
of physics, we may say that the true correlation of the physical 
■sciences lies in the equation of continuity 



9/ 






In the Theory of the Potential Function in the Attraction 
-of Gravitation, and in Electrostatics and Electrodynamics,* 
Fia Laplace's equation (3) is the value of the Potential Func- 
tion, at any external point {x^y, z\ due to any distribution of 
matter or of electricity; in the theory of the Conduction of 
Heat in a homogeneous solid f V is the- temperature at any 
point in the solid after the stationary temperatures have been 
■established, and in the theory of the irrotational flow of an 
incompressible fluid J F is the Velocity Potential Function 
and (3) is known as the equation of continuity. 

If we use spherical coordinates, (3) takes the form 
J. 9F\ 

* See Peirce*s Newtonian Potential Function. Boston, 
t See Fourier's Analytic Theory of Heat. London and New York, 1878 ; 
or Riemann's Partielle Differentialgleichungen. Brunswick. 

J See Lamb's Hydrodynamics. London and New York, 1895. 
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10 HARMONIC FUNCTIONS. 

and if we use cylindrical coordinates, the form 

In the theory of the Conduction of Heat in a homogene- 
ous solid,* u in equation (4) is the temperature of any point 
(;r, y, z) of the solid at any time /, and a^ is a constant deter- 
mined by experiment and depending on the conductivity and 
the thermal capacity of the solid. 

Art. 2. Homogeneous Linear Differential Equations* 

The general solution of a differential equation is the equa« 
tion expressing the most general relation between the primi- 
tive variables which is consistent with the given 'differential 
equation and which does not involve differentials or derivative^. 
A general solution will always contain arbitrary (i.e., undeter- 
mined) constants or arbitrary functions. 

A particular solution of a differential equation is a relation 
between the primitive variables which is consistent with the 
given differential equation, but which is less general than the 
general solution, although included in it. 

Theoretically, every particular solution can be obtained 
from the general solution by substituting in the general solu- 
tion particular values for the arbitrary constants or particular 
functions for the arbitrary functions ; but in practice it is often- 
easy to obtain particular solutions directly from the differential 
equation when it would be difficult or impossible to obtain the 
general solution. 

{a) If a problem requiring for its solution the solving of a 
differential equation is determinate, there must always begivea 
in addition to the differential equation enough outside condi- 
tions for the determination of all the arbitrary constants or 
arbitrary functions that enter into the general solution of the 
equation ; and in dealing with such a problem, if the differen- 
tial equation can be readily solved the natural method of pro- 
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HOMOGENEOUS LINEAR DIFFERENTIAL EQUATIONS. 11 

cedure is to obtain its general solution, and then to determine 
the constants or functions by the aid of the given conditions. 

It often happens, however, that the general solution of the 
dififerential equation in question cannot be obtained, and then, 
since the problem, if determinate, will be solved, if by any 
means a solution of the equation can be found which will also 
satisfy the given outside conditions, it is worth while to try to 
get particular solutions and so to combine them as to form a 
result which shall satisfy the given conditions without ceasing 
to satisfy the differential equation. 

{b) A differential equation is linear when it would be of the 
first degree if the dependent variable and all its derivatives 
were regarded as algebraic unknown quantities. If it is linear 
and contains no term which does not involve the dependent 
variable or one of its derivatives, it is said to be linear and 
homogeneous. 

All the differential equations given in Art. i are linear and 
homogeneous. 

{c) If a value of the dependent variable has been found 
-which satisfies a given homogeneous, linear, differential equa- 
tion, the product formed by multiplying this value by any 
-constant will also be a value of the dependent variable which 
will satisfy the equation. 

For if all the terms of the given equation are transposed 
to the first member, the substitution of the first-named value 
must reduce that member to zero ; substituting thp second 
value is equivalent to multiplying each term of the result of 
the first substitution by the same constant factor, which there- 
fore may be taken out as a factor of the whole first member. 
The remaining factor being zero, the product is zero and the 
equation is satisfied! ^ 

id) If several values of the dependent variable have been 
found each of which satisfies the given differential equation, 
their sum will satisfy the equation ; for if the sum of the values 
in question is substituted in the equation, each term of the sum 
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12 HARMONIC FUNCTIONS. 

will give rise to a set of terms which must be equal to zero, and 
therefore the suip of these sets mCist be zero. 

{e) It is generally possible to get by some simple device 
particular solutions of such differential equations as those we 
have collected in Art. i. The object of this chapter is to find 
methods of so combining these particular solutions as to satisfy 
any given conditions which are consistent with the nature of 
the problem in question. 

This often requires us to be able to develop any given func- 
tion of the variables which enter into the expression of these 
conditions in terms of normal forms suited to the problem with 
which we happen to be dealing, and suggested by the form of 
particular solution that we are able to obtain for the differential 
equation. 

These normal forms are frequently sines and cosines, but 
they are often much more complicated functions known as 
Legendre's Coefficients, or Zonal Harmonics ; Laplace's Coef- 
ficients, or Spherical Harmonics ; Bessel's Functions, or Cylin- 
drical Harmonics; Lamp's Functions, or Ellipsoidal Har- 
monics; etc. 

Art. 3. Problem in Trigonometrix: Series. 

As an illustration let us consider the following problem t 
A large iron plate tc centimeters thick is heated throughout 
to a uniform temperature of lOO degrees centigrade ; its faces 
are then suddenly cooled to the temperature zero and are kept 
at that temperature for 5 seconds. What will be the tempera- 
ture of a point in the middle of the plate at the end of that 
time? Given a^ =0.185 in C.G.S. units. 

Take the origin of coordinates in one face of the plate 
and the axis of X perpendicular to that face, and let u be the 
temperature of any point in the plate t seconds after the cooU 
ing begins. 

We shall suppose the flow of heat to be directly across the 
plate so that at any given time all points in any plane parallel 
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PROBLEM IN TRIGONOMETRIC SERIES. 13 

to the faces of the plate will have the same temperature. 

Then u depends upon a single space-coordinate x ; — = o and 

Pill 

— = o, and (4), Art. I, reduces to 

^ = a^?^. (I) 

Obviously, i u = ioo° when / = o, (2) 

u = o when x = o, (3) 

and u = o when x = tt; (4) 

and we need to find a solution of (i) which satisfies the con- 
ditions (2), (3), and (4). 

We shall begin by getting a particular solution of (i), and 
we shall use a device which always succeeds when the equa- 
tion is linear and homogeneous and has constant coefficients. 

Assume* u = ^^+>'^, where /3 and y are constants; substi- 
tute in (i) and divide through by ^^+'y' and we get y = a^/3^ ; 
and if this condition is satisfied, u = ^^*+y^ is a solution of (i). 

u = ^^+'''^'' is then a solution of (i) no matter what the 
value of /3. 

We can modify the form of this solution with advantage. 
Let /3 = /^/,t then u = ^-"VV*^* is a solution of (i\ as is also 

By (d), Art. 2, 

u = ^-^v A^ ^' 1 = ^-av«^ COS ;ix (5) 

is a solution, as is also 

^^^-a'^a/V f .^ = ^-«V^sin/i4^; . (6) 

and /4 is entirely arbitrary. 

* This assumption must be regarded as purely tentative. It must be tested 
by substituting in the equation, and is justified if it leads to a solution, 
t The letter % will be used to represent -f/"— i. 
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14 HARMONIC FUNCTIONS. 

By giving different values to /i we get different particular 
solutions of (i) ; let us try to so combine them as to satisfy our 
•conditions while continuing to satisfy equation (i). 

u = ^-*V*^ sin /ur is zero when ;r = o for all values of /4 ; it 
4s zero when ;r = ;r if /4 is a whole number. If, then, we write 
-u equal to a sum of terms of the form ^^ "*'*""' sin mx^ where 
m is an integer, we shall have a solution of (i) (see (d). Art. 2) 
^which satisfies (3) and (4). 

Let this solution be 

u =:A,e' -•* sin x + A,e '^^' sin 2x + A^e ' »*'* sin 3:1: + ... , (7) ^ 

vij, A^, A^y . . . being undetermined constants. 
When / = o, (7) reduces to 

u = A^ sin X -}- A^ sin>2ar + -^1 sin 3:1: + • • • • (8) 

If now it is possible to develop unity into a series of the 
form (8) we have only to substitute the coefficients of that 
series each multiplied by 100 for A^, A^, ^,. . . in (7) to have 
a solution satisfying (i) and all the equations of condition (2), 
<3), and' (4). 

We shall prove later (see Art. 6) that 

I = — sin ;ir + " sin 3;r + 7 sin 5;i: + • • • 
^L 3 5 -J 

'for all values of x between o and n. Hence our solution is 
u = ^r^""'' sin X + -e'^'' sin ^x + ^e"^'' sin S^r + . . . | (9) 

To get the answer of the numerical problem we have only 

to compute the value of u when ;r = — and / = 5 seconds. As 

there is no object in going beyond tenths of a degree, four- 
place tables will more than suffice, and no term of (9) beyond . 

the first will affect the result. Since sin — = i, we have to 

•compute the tirumerical value of 

• Digitized by VjOOQIC 



PROBLEM IN ZONAL HARMONICS. 15 

400 

^—e'"'^^ where a^ = 0.185 and / = 5. 

log a^ = 9.2672 — 10 log 400 = 2.6021 

log t = 0.6990 colog n = 9.5059 — 10 

log €^t = 9.9662 — 10 colog ^*^ = 9.5982 — 10 

log log. =9,6378-10 ^-^^ 

log log ^'' = 9.6040 — 10 

log ^'^ = 0.4018 u = 5o''«8 . 

If the breadth of the plate had been c centimeters instead 
of n centimeters it is easy to see that we should have needed 
the development of unity in a series of the form 

A^ sm V ^% sm 1- A^ sin + 

Prob. I. An iron slab 50 centimeters thick is heated to the tem- 
perature 100 degrees Centigrade throughout. The faces are then sud- 
denly cooled to zero degrees, and are kept at that temperature for 
10 minutes. Find the temperature of a point in the middle of the 
slab, and of a point 10 centimeters from a face at the end of that 
time. Assume that 

4/ . 7TX , I . XTtX , I . KTtx , \ ^ 
I = - sin h " sin h - sin \-.. . I from x = o to x = c. 

Ans, 84^.0; 49°.4. 

Art. 4. Problem in Zonal Harmonics. 

As a second example let us consider the following problem : 
Two equal thin hemispherical shells of radius unity placed 
together to form a spherical surface are separated by a thin 
layer of air. A charge of statical electricity is placed upon 
one hemisphere and the other hemisphere is connected with 
the ground, the first hemisphere is then found to be at poten- 
tial I, the other hemisphere being of course at potential zero. 
At what potential is any point in the " field of force" due to 
the charge? 

We shall use spherical coordinates and shall let Fbe the 
potential required. Then Fmust satisfy equation (5), Art. i. 
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16 HARMONIC FUNCTIONS. 

But since from the symmetry of the problem V is obviously 
independent of 0, if we take the diameter perpendicular to the 

plane separating the two conductors as our polar axis, -— j is 
zero, and our equation reduces to 

ar" "T" sin d d^ ' 

Fis given on the surface of our sphere, hence 

V=:f{ff) when r=i, (2) 

where y(6^) = I if o < » < -, and /{ff) = oii-<0<7r. 

2 2 

Equation (2). and the implied conditions that V is zero at 
an infinite distance and is nowhere infinite are our conditions. 

To find particular solutions of (i) we shall use a method 
which is generally efifective. Assume* that V=R@ where R 
is a function of r but not of 6, and is a function of ff but 
not of r. Substitute in (i) and reduce, and we get 

2 rd\rR) ^ i_ 4"^" ^-dBJ . (3) 

R dr" & sin 6 dO 

Since the first member of (3) does not contain d and the 
second does not contain r and the two members are identically 
equal, each must be equal to a constant. Let us call this 
constant, which is wholly undetermined, w(w+ i) ; then 

J d&\ 

R-di^ = -&^ire — Te «'('«+ 0. 

d'irR) 
whence r \ , — m(m -\-\)R = o, (4) 

^"*^ -^re—de — + '«('« +i)® = o. (5) 

* See the first foot-note on page 175. 
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PROBLEM IN ZONAL HARMONICS. 17 

Equation (4) can be expanded into 

^d-R ^ dR , , ,^ 

r«— + 2r-^ - m(fn-\-\)R = o, 

and can be solved by elementary methods. Its complete 
solution is 

^ = ^r'« + ^r-«-^ (6) 

Equation (5) can be simplified by changing the independ- 
ent variable to x where x = cos 6, It becomes 

dV d@l 

an equation which has been much studied and which is known 
as Legendre's Equation. 

We shall restrict w, which is wholly undetermined, to posi- 
tive whole values, and we can then get particular solutions of 
(7) by the following device : 

Assume* that can be expressed as a sum or a series of 
terms involving whole powers of x multiplied by constant 
coefficients. 

Let = ^anX"" and substitute in (7). We get 

^\n{n — i)^«;ir*-' — n{n + i>„jr« + m{m + i)^.^*] = o, (8) 

where the symbol ^ indicates that we are to form all the 
terms we can by taking successive whole numbers for n. 

Since (8) must be true no matter what the value of x, the 
coefficient of any given power of x, as for instance x^y must 
vanish. Hence 

{k + 2Xk + i)^;k+, - k{k + i)^^ + m{m + i)af, = o, 
and «».. = - (^_^,)(^_^,| ' a,. (9). 

If now any set of coefficients satisfying the relation (9) be 
taken, = ^Uj^ will be a solution of (7). 

If k = m, then Uj,^^ = o, ^a;*+, = o, etc. 

* See the first foot-note on page 175. 
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18 HARMONIC FUNCTIONS. 

Since it will answer our purpose if we pick out the simplest 
set of coefficients that will obey the condition (9), we can take 
a set including a^. 

Let us rewrite (9) in the form 

«* - - ^jn-k){m + k-x)' ^^°> 

, We get from (lo), beginning with k=.m — 2, 

m{m — i) 

^ _ m{m — i)(w - 2){m — 3) ^ 
*"'* 2.4.(2^— 1X2W— 3) *"' 

_ _ m{m ~ i){m - 2){m - 3){m - 4){m - 5) 
'^"— " 2 . 4 . 6 . (2w - i){2m - 3X2W - 5) '''" ' ^^^* 

If m is even we see that the set will end with a^; if m 
is odd, with a,. 

*" L 2. {2m— I) 

"^ 2.4.(2W— l)(2W— 3) '"J 

where ^3!« is entirely arbitrary, is, then, a solution of (7). It is 
found convenient to take a^ equal to 

(2m — i)(2;« — 3) ... I 

and it will be shown later that with this value of ^j,„ , — i 
when X = I, 

is a function of x and contains no higher powers of x 
-than x"^. It is usual to write it as -P^C^r). 

We proceed to write out a few values of Pm{^) from, the 
formula 

p (^\ = (2w - i)(2^ - 3) ' - I r^ _ m(m- i) ^^., 
'"^ ^ ;«! L 2.(2w— i) 

~ 2.4.(2W- 1X2W — 3) J 
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(12) 



We have : 

Plx) = I or P.(cos 0) = i, 

Pi(x) = X or PJ^cos ff) = cos 0, 

P,{^) = 4(3^* - I) or ^,(cos ff) = Us cos'</ - i), 

P,{x) = Hs^' - 3x) or /'.(cos 0) = i(5 cos'6' - 3 cos 0), 

Pix) = i(35jr* - sox' + 3) or 

/'/cos 0) = K35 cos*^' - 30 cos'6' + 3), 
^.W = i(63^* - 70X' + 15^:) or 

/'.(cos 6/) = ^(63 cos' — 70 cos'^+ 15 cos 6*). , 

We have obtained © = /'„(jr) as a particular solution of (7), 
and © = P„(cos ^ as a particular solution of (5). Pj{x) or 
/'^(cos 61) is a new function, known as a Legendre's Coefficient, 
or as a Surface Zonal Harmonic, and occurs as a normal form 
in many important problems. 

V= r''P^{cos 0) is a particular solution of (i), and r*'/'„(cos tf) 
is sometimes called a Solid Zonal Harmonic. 

V = A,P,{cos 0) + A.rPScos 0) + A^PJ^cos 0) 

+ AyP,{cos0)-\-... (13) 
satisfies (i), is not infinite at any point within the sphere, and 
reduces to 
V= A,P,{cos 6/) + ^,P.(cos 0) + A,P,{cos 0) 

+ ^./'.(cos6l)+. 



when r = I. 

r^_ ^.^.(cosg) , AA(cos0) A,P,(cos0) 

A,P,icos 0) 



+ 



(14) 



(15) 



satisfies (i), is not infinite at any point without the sphere, is 
equal to zero when r = oo , and reduces to (14) when r = i. 

If then we can develop /(^ [see eq. {2)] into a series of the 
form (14), we have only to put the coefficients of this series in 
place of the A^, A^, A^y ... in (13) to get the value of Ffor a 
point within the sphere, and in (15) to get the value of Fat a 
point without the sphere. 
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We shall see later (Art. i6, Prob. 22) that if /{O) = i for 
iO < e<—- 2Xidf{0) = ofor—-<e<7t, 

2 2 



+77'2-i|^'(~^'>-- ^'^ 



Hence our required solution is 



"r'-^—Zncosff)^-^^ (17) 



at an internal point ; and 

.3 
4 



^=i+!i^^«(~^^-r^-^^'(^^^^ 



. li.illip^ 3^_ _ (i8) 

' 12 2.4r' •^ ^ ^ ' 



?t an external point. 

If r = — and 6^ = o, (17) reduces to 

F=-+~ 5- tH — ^- -T-..., sncei^«(i) = i. 

2*44 8 2 4" ' 12 2.4 4» ' «\ / 

To two decimal places V = 0.68, and the point r = -, 6^ = o 

4 
is at potential a68. 

If r = 5 and 6^ = — , (18) and Table I, at the end of this 
4 

chapter, give 

^=^+!'T'-°'7°7i + |- ^•-^•0.1768+... =0.12, 

and the point r = $y d •= is at potential 0.12. 

4 

If the radius of the conductor is a instead of unity, we 

iiave only to replace r by - in (17) and (18). 

a 
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PROBLEM IN BESSEL's FUNCTIONS. '4,1 

Prob. 2. One half the surface of a solid sphere 12 inches in di- 
ameter is kept at the temperature zero and the other half at 100 de- 
grees centigrade until there is no longer any change of temperature 
at any point within the sphere. Required the temperature of the 
center ; of any point in the diametral plane separating the hot and 
cold hemispheres ; of points 2 inches from the center and in the 
axis of symmetry ; and of points 3 inches from the center in a di- 
ameter inclined at an angle of 45° to the axis of symmetry. 

Ans. 50°; so'*; 73^.9; 26°.! ; 77^1; 22^9. 

Art. 5. Problem in'Bessel's Functions. 

As a last example we shall take the following problem : 
The base and convex surface of a cylinder 2 feet in diameter 
and 2 feet high are kept at the temperature zero, and the upper 
base at 100 degrees centigrade. Find the temperature of a 
point in the axis one foot from the base, and of a point 6 inches 
from the axis and one foot from the base, after the permanent 
state of temperatures has been set up. 

If we use cylindrical coordinates and take the origin in the 
base we shall have to solve equation (6), Art. i ; or", represent- 
ing the temperature by u and observing that from the sym- 
metry of the problem u is independent of 0, 

subject to the conditions 

u = when js = Oy (2) 

• « = o " r= I, (3) 

U = 100 . " £r = 2. (4) 

Assume u = RZ where R is a, function of r only and Z of 
js only; substitute in (i) and reduce. 

„, , I {PR , I dR I (PZ 

The first member of (5) does not contain z ; therefore the 
second member cannot. The second member of (5) does hot 
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contain r; therefore the first member cannot. Hence each 
member of (5) is a constant, and we can write (5) 

Rdr^'^rRdr" Z ds^ ^ ^' ^^ 

when /I* is entirely undetermined. 

Hence ^-;,«Z=o, (7) 

Equation (7) is easily solved, and its general solution is 

Z = Ae**-* '\- Be'^, or the equivalent form 

Z-=- C cosh {^}x£) + Z? sinh {ptz). (9) 

We can reduce (8) slightly by letting //r = x^ and it becomes 

d^R , I dR . r, / N 

Assume, as in Art. 4, that /? can be expressed in terms of 
whole powers of x. Let R = 2a^* and substitute in (10). 
We get 

2[n{n - i)a^x- ' • + «^«;r«- ' + a^^] = o, 

an equation which must be true, no matter what the value of x. 
The coefficient of any given power of ;r, as;r*"', must, then, 
vanish, and 

^(^ — iK + kaj^ -^aj^,, = o, 
or ^*^* + ^*-, = o, 

whence we obtain a^^., = — i^a^ (11) 

as the only relation that need be satisfied by the coefficients in 
order that R = 2aj^ shall be a solution of (10). 

If ^ = o, ^i_ , = o, a^.^ = o, etc. 

We can, then, begin with ^ = o as the lowest subscript. 
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From (if) a*= — 



e 



Then «. = -^. «.= ^. ''• = -?:7:6"*'^- 

Hence ie = a.[i - f; + -^. - -j-^. + . .. ], 

where a, may be taken at pleasure, is a solution of (lo), pro- 
vided the series is convergent. 

Take a^^\, and then R z=.JJ^x) where 

/oW = ^ - 1.+^. - ?:^«+2'.4'.6".y "■ • ' • ^^^^ 

is a solution of (lo). 

Jlx) is easily shown to be convergent for all values real or 
imaginary of x, it is a new and important form, and is called a 
BesseFs Function of the zero order, or a Cylindrical Har- 
monic. 

Equation (lo) was obtained from (8) by the substitution of 
X := pir\ therefore 

is a solution of (8), no matter what the value of ;i; and 
u = ///ir) sinh (/ta) and u = /,(//r) cosh (//^r) are solutions of 
(i). u = JS.fJir) sinh (//^) satisfies condition (2) whatever the 
value of //. In order that it should satisfy condition (3) /i 
must be so taken that 

Uti) = 0', (13) 

that is, /i must be a root of the transcendental equation {13). 

It was shown by Fourier that Jlji) = o has an infinite num- 
ber of real positive roots, any one of which can be obtained to 
any required degree of approximation without serious diffi- 
culty. Let /i, , //, , /i, , . . . be these roots ; then 

u = AJl^y) sinh (jx,z) + AJl^i^r) sinh {^^) 

+ AJlfi,r) sinh (/i,^) + . . . (14) 

is a solution of (i) which satisfies (2) and (3). 
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24 HARMONIC FUNCTIONS. 

If now we can develop unity into a series of the form 
I = BJXf^^r) + BJJ^/) + ^./.(/i.r) + . . . , 

satisfies (i) and the conditions (2), (3), and (4). 

We shall see later (Art. 21) that \i /.{x) = - ^^ 

ax 

for values of r < i. 
Hence 

Ly«./.(A'.) sinh (2//,) /',/.(/',) smh (s/'O -J 

is our required solution. 

At the point r = 0, z= i (17) reduces to 

u = ^oof ^'"^^' + —E}hJh |_ 1 

L/^./.(a<,) sinh (2//,) //,/,0^,) sinh (2/*,) ^ 'J 



= 100 - 



+ 1 + 1 

/^./,(/^.) cosh /i, J' 



L/^i/iO^i) cosh /i, /^,/,(/^,) cosh /i, 
since /^(o) = i and sinh {2x) = 2 sinh ;r cosh x. 

If we use a table of Hyperbolic functions* and Tables II 
and III, at the end of this chapter, the computation of the 
value of u is easy. We have 

/i, = 2.405 //,= 5.520 

/i(/^i) = 0.5 190 /,{m,) = - 0.3402 



colog /ij = 9.6189 — 10 colog /ij = 9.2581 — 10 
- /Xa) = 0.2848 " /U)= o.4683« 

" cosh /!,= 9.2530 — ID " cosh/i,= 7.9037 — 10 

9.1567 — 10 • 7.6301 « — 10 

* See Chapter IV, pp. 162, 163, for a four- place table on hyperbolic func- 
tions. 
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PROBI EM IN BESSEL's FUNCTIONSw 25 

0"i/xO^i) cosh /I,)'* = o. 1434 
(/^/i(^«) cosh /I,)-* = - 0.0058 

0.1376; u=if.S 

At the point r = J, ^ = i, (17), reduces to 

L A'./iC/^,) cosh Ml M,/i{M,) cosh /i, ' ' J* 
/,(i/i.) = 0.6698 

log/.(i/iO= 9-8259 -10 
colog /ij/,(/i,) cosh /i, = 9.1 567 — 10 

8.9826 — 10; 
/o(i/^.) = ~ 0.1678 



log AHm.) = 9-2248« - 10 
colog Mt/iiMi) cosh /i, = ' 7.630i« — ID 

6.8549 — 10; 



/^,/i(/^i) cosh/i, 



= 0.0961 



/oiJM.) ^ 0.0007 , 

;^«/,{yw«) cosh M, 0.0968 ' « = 9^7 

If the radius of the cylinder is a and the altitude ^, we have 
only to replace jdhy pia in (13) ; 2/^^, 2//,, ... in the denonni- 
nators of (15) and (17) by //j^, //,^, . . . ; and pi^^ //,, ^u,, ... in 
the denominators of (16) and (17) by /i^a, pi^a, jn^a, .... 

Prob. 3. One base and the convex surface of a cylinder 20 cen- 
timeters in diameter and 30 centimeters high are kept at zero tem- 
perature and the other base at 100 degrees Centigrade. Find the 
temperature of a point in the axis and 20 centimeters from the cold 
base, and of a point 5 centimeters from the axis and 20 centimeters 
from the cold base after the temperatures have ceased to change. 

Ans. 1 3°. 9; 9°.6. 
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Art. 6. The Sine Series. 

As we have seen in Art. 3, it is sometimes important to be 
able to express a given function of a variable, x^ in terms of sines 
of multiples of x. The problem in its general form was first 
solved by Fourier in his ** Th^orie Analytique de la Chaleur** 
(1822), and its solution plays an important part in most branches 
of Mathematical Physics. 

Let lis endeavor to so develop a given function of x,f{x\ 
in terms of sin x, sin 2;r, sin ix, etc., that the function and the 
series shall be equal for all values of x between o and n. 

We can of course determine the coefficients ^,, ^,, ^„ . . . a^ 
so that the equation 

f{pc) = a^ sin X '\' a^ sin 2jr + ^« sin jjr 4- • • • + ^« sin nx (i) 
shall hold good for any n arbitrarily chosen values of x between 
O and n ; for we have only to substitute those values in turn) 
in (i) to get n equations of the first degree, in which the n co- 
efficients are the only unknown quantities. 

For instance, we can take the n equidistant values Ax^ 2Axy. 

%Ax, , . . tiAx. where Ax = — ; — , and substitute them for;*: ia 

(i). We get 

f{Ax) = a, sin Ax + ^, sin 2Ax + ^, sin ^Ax + . . . 

+ a« sin nAXf 

f{2Ax) = a^ sin 2Ax + ^, sin /^Ax + a^ sin 6 Ax + • • • 

+ a^sxn 2nAx^ 

fij^Ax") = a^ sin ^Ax + ^, sin 6Ax + a^ sin ^Ax + . . . 

+ a^ sin ytAxy 

/{nAx) = a^ sin nAx + ^, sin 2nAx + a, sin ^nAx + . . . 

+ an sin n*Axy 
n equations of the first degree, to determine the n coefficients 

^j » ^, » ^f > • • • ^w 

Not only can equations (2) be solved in theory, but they 
can be actually solved in any given case by a very simple and 



(2> 
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THE SINE SERIES. 27 

ingenious method due to Lagrange,* and any coefficient a^n can 
be expressed in the form 

2 *"* 

a^ = Z-TTj^ A^^^) sin {Km^x). (3) 

If now n is indefinitely increased the values of x for which 
{1) holds good will come nearer and nearer to forming a con- 
tinuous set ; and the limiting value approached by a^ will 
probably be the corresponding coefficient in the series required 
to represent /{x) for all values of x between zero and zr. 

Remembering that (« + 0-^-^ = ^> the limiting value in 
question is easily seen to be 

V 

a^ = -CA^) sin mxdx. (4) 


This value can be obtained from equations (2) by the fol- 
lowing device without first solving the equations : 

Let us multiply each equation in (2) by the product of /4x 
and the coefficient of a^ in the equation in question, add the 
equations, and find the limiting form of the resulting equation 
as n increases indefinitely. 

The coefficient of any a, a^ in the resulting equation is 

sin kAx sin mAx. Ax -{- sin 2k Ax sin 2m Ax .Ax-\-,.. 
+ sin nKAx sin nmAx. Ax. 

Its limiting value, since (« + i)Ax = a-, is 



IT 

/ sin Kx sin mx.dx; 



but 



f sin KX sin mx . dx =i\J\(:os {m — k)x — cos{m + K)x'\dx=^o 

« 

if m and k are not equal. 

* See Riemann's Pariielle Diflfercntialglcichungen, or Byerly's Fourier's 
Series and Spherical Harmonics. 
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The coefficient of a^^ is 
J;i;(sin* wJjr+sin" 2wJ;t: + sin' 3wJ;r+ . . . +sin' nmAxy 
Its limiting value is 

/ sin' mx .dx = -. 



I'he first member is 

f{^x) sin m^x . /ix -\-/{2^x) sin 2inAx . J;r + • • • 
■\-f(nAx^ sin mnAx . J;r, 
and its limiting value is 

/ y(;ir) sin mx . ^;ir. 



Hence the limiting form approached by the final equation 
as n is increased is 



/ f{x) sin mx . dx = —a^. 



Whence a^ = -^ A^) sin mx.dx (5) 

as before. 

This method is practically the same as multiplying the 
equation 

f{pc) = a^ sin ^ + ^j sin 2;r + ^i:, sin 3;ir + • . . (6) 

by sin mx >dx and integrating both members from zero to n. 

It is important to realize that the considerations given in 
this article are in no sense a demonstration, but merely estab- 
lish a probability. 

An elaborate investigation * into the validity of the develop- 
ment, for which we have not space, entirely confirms the results 
formulated above, provided that between ;r = o and x^=z n the 

* See Art. 10 for a discussion of this question. 
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function is finite and single-valued, and has not an infinite num- 
ber of discontinuities or of maxima or minima. 

It is to be noted that the curve represented by y:=.f{x) 
need not follow the same mathematical law throughout its 
length, but may be made up of portions of entirely different 
curves. For example, a broken line or a locus consisting of 
finite parts of several different and disconnected straight lines 
can be represented perfectly well by ^ = a sine series. 

As an example of the application of formula (5) let us take 
the development of unity. 

Here f{x) =1. 



IT 

c^m^ - j sin mx . dx ; 



/ 



, cos mx 

sin ntx .dxz=z . 

m 



It 



sin mx.dx = — (i — cos mii) = — fi — (— i)*"! 

m^ ' m ^ ' ^ 

= o if /« is even 

= — if w is odd. 
m 



Hence I =l(^ + ?i^ii£ + ^+?iil^ + ...). (7) 

It is to be noticed that (7) gives at once a sine development 
for any constant c. It is, 

^^^^sin_£^sin35^sin5£_^ j (g) 

Prob. 4. Show that for values of x between zero and n 

i X Fsin x sin 2x , sin xx sin 4^; , "1 

(.) .= ,L_-_- + _3___4f + ...J, 

lK\ f(^\ — 4 r?!!Lf sinj* sin 5a: sin ^x "I 

W/(*) - -\-pr - -Y' + ~' f- + '"] 



Digitized by 



Google 



30 HARMONIC FUNCTIONS, 

n It 

\if(x) =■ xioxo <x < -^ and f{pc) -^^ n — x iox - < x < n. 

w /(*) = 

2 Fsin X 2 sin 2X , sin 3^ , sin 5^ , 2 sin 6jf , sin 7^ , "1 

if /(^) = I for o < :c < -, and f{x) = o for — < ^ < «". 

(^) sinh :c = 

2 sinh zrri. 2. i3- 4- , "1 

— sin ^ sin 2^ + sm 3^ sin 4^ + . . . . 

zr L2 5 ' 10 '^ 17 ^ ^ J 

ilvT" " ?) ^'"^ "^ "" T ^'"^ ^^ "^ (l" " 3"")^^'' ^^"^ ^""'^ 4^+ • • •]• 

Art. 7. The Cosine Series. 

Let us now try to develop a given function of at in a series 
of cosines, using the method suggested by the last article. 
Assume 

f{x) = ^^ + ^, cos X'\-b^ cos 2x + b^ cos ix-\- . , . (i) 

To determine any coefficient b^ multiply (i) by cos mx .dx 
.and integrate each term from o to n. 

ir 

jb^ cos mx .dx = 0, 



w 

/ bk cos kx cos w;ir . ^;r =0, if w and ^ are not equal. 


w 

/ b^ cos" »^;t: .^;ir = — ^„, if w is not zero. 



Hence b^^ — J /{x) cos mx . dx, (2) 


if m is not zero. 
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To get'^, multiply (i) by dx and integrate from zero to n. 
Jb,dx = b^n, 



/ bk cos kx.dx = 0. 



rr 

Hence b,= 'Lff{x)dx, (3) 



which is just half the value that would be given by formula (2) 
if zero were substituted for m. 

To save a separate formula (l) is usually written 

f{x) = i^, + b^ cos x + b^ cos 2x + b^ cos 3^ + • • •» (4) 

and then the formula (2) will give b^ as well as the other coef- 
ficients. 

Prob. 5. Show that for values of x between o and 7t 

, V 7t A /cos X , cos ^x . cos 5^^ , \ 

if f(x) := xiox o < X < -, and f{x) =^ 7t — x for - < x X ^; 

2 2 

/ . r/ \ I I 2 /cos ^ cos 3a: , cos 5:*: \ 

(^)A-)=-+^[- s+-r -■■•)' 

It 7t 

if f{pc) = I for o < a: < -, and f{pc) =^ oiox — < x <n^ 

2 2 

2 I I I 
(d) sinh ^ = - —(cosh tt — i) (cosh ^r + i) cos ^ 

+ -(cosh 7t — i) cos 2A: (cosh ;r -|- i) cos 3^: + . . . ; 

, V a 7t^ /cos X cos 2^ . COS %X COS 4:V . \ 
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Art. 8. Fourier's Series. 

Since a sine series is an odd function of x the development 
of an odd function of x in such a series must hold good from 
X •=> — n \.o X •= n^ except perhaps for the value ^ = o, where 
it is easily seen that the series is necessarily zero, no matter 
what the value of the function. In like manner we see that 
if J\x) is an even function of x its development in a cosine 
series must be valid from ;r = — ;r to ;r = ;r. 

Any function of x can be developed into a Trigonometric 
series to which it is equal for all values of x between — n and n. 

\^t\.f{x) be the given function of x. It can be expressed 
as the sum of an even function of x and an odd function of x 
by the following device : 

^^) ^ A^) +/(- ^) I A^) -A- X) (,) 

identically; but ^ ^ ' "^ is not changed by reversing 

the sign of x and is therefore an even function of x\ and when 

we reverse the sign of ;r, -^—^ — --^ is affected only to the 

extent of having its sign reversed, and is consequently an odd 
function of x. 

Therefore for all values of x between — n and n 

A^) +/(- ^) ^ i^^ _|_ ^^ cos^r + b^ cos 2x + b, cos 3;r + . . . 

2 2 

where b^= — I '- cos mx.dx\ 

ltd 2 



and '^-^^-^ — ^^ = ^1 sm :r + a, sm 2x + a^ sm j;r + • . . 

2 

2 /» fix\ /T x\ 

where ^« = - / "^^^-^ — — sin mx . dx. 

nJ 2 
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b^ and a^ can be simplified a little. 

7t J 2 



W IT 

= — y/W COS mx . dx'\-jj\-^x) cos /«;r . ^;r |; 
but if we replace ;t: by — ;r, we get 
J A— ^) COS mx.dx^i^J f[x) cos mx.dxz=z J /[x)cos mx.dx, 

_» 

and we have ^^ = — j f(x) cos wjr . ^/;r. 

In the same way we can reduce the value of a^ to 

— I f[x\ sin mx .dx. 
Hence 

/(;ir) = - ^0 + ^1 cos ;ir + ^t cos 2^: + ^. cos 3;r + • • • 

+ ^, sin ^ + a;, sin 2x -{- a^sin 3^ + . . . , (2) 

where *« = —jA^) cos w;ir . dx, (3) 



and 



^«, = -^JA^) sin ^-^ • dx, (4) 



and this development holds for all values of x between — n 
and n. 

The second member of (2) is known as a Fourier's Series. 

The developments of Arts. 3 and 7 are special cases of 
development in Fourier's Series. 

Prob. 6. Show that for all values of x from — n \o n 

2 sinh TrFi i ,1 i ,1 , "I 

n \_2 2 5 ip '^ 17 ' I 
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, 2 sinh TcFi , 2 . , 3 • 4 • , "1 

H — sin jc sin 2JIP + -^ sin 3^ sin 4;c + . . . . 

^ L2 5 10 '^ 17 J 

Prob. 7. Show that formula (2), Art. 8, can be written 

/{x) = i C^ cos/?, + C, cos {X - /? J + C^ cos {2X - /?,) 

+ ^, cos (3^ - A) + . • • T 
where ^« = (««.* + ^..")* and /?« = tan'^ |^- 

Prob. 8. Show that formula (2), Art. 8, can be written 

/{x) =^c, sin /?o + c, sin (^ + /? J + ^a sin (2^ + /?,) 
2 

+ ^3 sin (3^ + /?J + . . . , 
where ^« = K' + O* and /?« = tan'^^. 



Art. 9. Extension of Fourier's Series. 

In developing a function of x into a Trigonometric Series it 
IS often inconvenient to be held within the narrow boundaries 
X = — 7C and ;r = tt. Let us see if we cannot widen them. 

Let it be required to develop a function of x into a 
Trigonometric Series which shall be equal to /{x) for all values 
of X between x = — c and x = c. 

Introduce a new variable 

2 = — X. 
c 

which is equal to — ;r when ;r = — ^, and to n when x =^ c. 

f{x) = /( — 2] can be developed in terms of js by Art. 8, 
(2), (3), and (4). We have 

+ a^ sin z -\-a^ sin 2^: + ^3 sin 3^ + . . . , (i) 
where b^ = — I f\ — ^] ^^^ ^^^ • ^^^ (2) 
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and a^ = —ff\^^ sin mz . dz, (3) 

— » 

and (i) holds good from ^ = — tt to ^ = ;r. 

Replace z by its value in terms of x and (i) becomes 

+ ^, sm-^+^,sm-^ + ^,sm-y-+...; (4) 
and (2) and (3) can be transformed into 

i.n = -yffi^) cos.^^dx, (5) 

a„ = -LJ'/{x) sin ^^dx, (6) 

—c 

and (4) holds good from ;ir = — ^ to ;r = ^. 

In the formulas just obtained c may have as great a value 
as we please so that we can obtain a Trigonometric Series for 
/{x) that will be equal to the given function through as great 
an interval as we may choose to take. 

It can be shown that if this interval c is increased indefi- 
nitely the series will approach as its limiting form the double 

00 00 

integral — I f{X)dX I cos a{\ — x)da^ which is known as a 

—00 

Fourier's Integral. So that 

+00 ao 

f{x) = -^f f{X)d\ f cos a{X - x)da (7) 

-09 . 

for all values of x. 

For the treatment of Fourier's Integral and for examples 
of its use in Mathematical Physics the student is referred to 
Riemann's Partielle Dififerentialgleichungen, to Schlomilch's 
Hohere Analysis, and to Byerly's Fourier's Series and 
Spherical Harmonics. 
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Prob. 9. Show that formula (4), Art. 9, can be written 
/(*) = ^^. cos /?, + c, cos (7^ - A ) + ^, cos \^^ - /?,) 

+ ..cos(2?^-A)+..., 

where c, - (^.; ^ ^«')* and /S^ = tan"^ -~^. 

Prob. 10. Show that formula (4), Art. 9, can be written 
/(^•) = \c, sin /?, + c, sin (^ + /? J + r, sin ^^ + /?,j 

+ ..sin(^ + A)+..., 
where ^« = (dr«' + O* and /?« = tan-i-^. 

Art. 10. Dirichlet's Conditions. 

In determining the coefficients of the Fourier's Series rep- 
resenting f{x) we have virtually assumed, first, that a series of 
the required form and equal to f(x) exists ; and second, that 
\\, \s uniformly convergent ; and consequently we must regard 
the results obtained as only provisionally established. 

It is, however, possible to prove rigorously that if /(^) is 
finite and single-valued from ;r=— ;rto;r=;r and has not 
an infinite number of (finite) discontinuities, or of maxima or 
minima between x ^= — n and x -=• n^ the Fourier's Series of 
(2), Art. 8, and that Fourier's Series only, is equal to f{x) 
for all values of x between — n and n, excepting the values of 
X corresponding to the discontinuities of /(;r), and the values 
n and — n ; and that if c is a value of x corresponding to a 
discontinuity of f(x\ the value of the series when ;r = ^ is 

i^J^*J^[y(^+e)+/(^-€)]; and that when ;r = ;r or 
X ^ —n the value of the series is ^[/(^r) +/(— ;r)]. 

This proof was first given by Dirichlet in 1829, and may be 
found in readable form in Riemann's Partielle Differential- 
gleichungen and in Picard's Traits d*Analyse, Vol. I. 
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A good deal of light is thrown on the peculiarities of trigo- 
nometric series by the attempt to construct approximately the 
curves corresponding to them. 

If we construct y ^ a^ sin x and j = ^, sin 2x and add the 
ordinates of the points having the same abscissas, we shall ob- 
tain points on the curve 

^ = ^j sin ;ir + a^ sin 2x, 

If now we construct ^ = ^, sin ^x and add the ordinates to 
those of J = ^, sin x -{■ a^ sin 2x we shall get the curve 

y =. a^sxw X -{- a^ sin 2x + a^ sin 3;r. 

By continuing this process we get successive approximations to 

y ■=. a^sm X ■\- a^ sin 2x + a^ sin 3jr + ^^ sin 4:1: + • • • 
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Let us apply this method to the series 
J = sin ;r + I sin 3;ir + I sin Sjt + (i) (See (7), Art. 6.) 

y = o when jr = o, — from ;r = o to jr = ;r, and o when x = 7t. 
4 

It must be borne in mind that our curve is periodic, hav- 
ing the period 2;r, and is symmetrical with respect to the 
origin. 

The preceding figures represent the first four approxima- 
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tion to this curve. In each figure the curve y = the series, 
and the approximations in question are drawn in continuous 
lines, and the preceding approximation and the curve corre- 
sponding to the term to be added are drawn in dotted lines. 

Prob. II. Construct successive approximations to the series 
given in the examples at the end of Art. 6. 

Prob. 1 2. Construct successive approximations to the Maclaurin's 

Series for sinh x^ namely x -{ — p H j- H 

Art, 11. Applications of Trigonometric Series. 

{a) Three edges of a rectangular plate of tinfoil are kept 
at potential zero, and the fourth at potential I. At what po- 
tential is any point in the plate ? 

Here we have to solve Laplace's Equation (3), Art. i, 
which, since the problem is two-dimensional, reduces to 

subject to the conditions V = o when ;r = o, (2) 

V=o ^^ X = a, (3) 

F = o " J = o, (4) 

F= I " y = d. (5) 

Working as in Art. 3, we readily get sinh /3y sin fijv, 
sinh)6[;/ cos/3x, cosh/fifj sin /3Xy and coshy6[;/ cos /3x as particu- 
lar values of F satisfying (i), 

V = sinh — ^ sin -^ satisfies (i), (2), (3), and (4). 

psinh ^ sinhi^ n 

V = M ^sin^+l ^,sin3JL^+... (6) 



is the required solution, for it rediices to i when y ^=. b. See 
(7), Art. 6. 
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(b) A harp-string is initially distorted into a given plane 
curve and then released ; find its motion. 

The differential equation for the snnall transverse vibrations 
of a stretched elastic string is 

as stated in Art. i. Our conditions if we take one end of 
the string as origin are 

(2) 

(3) 

(4) 

(5) 

Using the method of Art. 3, we easily get as particular solutions 
of (I) 

y = sin /3x sin a/3^, y = sin px cos apt, 

y = cos fix sin a fit, and y = cos fix cos a fit. 

. mndc mnat . ^ 
y = sin —J- cos -J- satisfies (i), (2), (3), and (4)1 

IJir ■ . mux mnat tfx 

y = ^a^ sm —y— cos — ^— , (o> 



where a„ = jCA^) sin ^j^ • ^^ (7> 



is our required solution ; for it reduces to f{x) when / = o. See 

Art. 9. 

Prob. 13. Three edges of a square sheet of tinfoil are kept at 
potential zero, and the fourth at potential unity ; at what potential 
is the centre of the sheet ? Ans. 0.25. 

Prob. 14. Two opposite edges of a square sheet of tinfoil are 
kept at potential zero, and the other two at potential unity ; at 
what potential is the centre of the sheet ? Ans. 0.5. 

Prob. 15. Two adjacent edges of a square sheet of tinfoil are 
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.kept at potential zero, and the other two at potential unity. At 
what potential is the centre of the sheet ? Ans. 0.5. 

Prob. 16. Show that if a point whose distance from the end of a 

harp-string is -th the length of the string is drawn aside by the 
n 

player's finger to a distance b from its position of equilibrium and 

then released, the form of the vibrating string at any instant is given 

.by the equation 

2bt^ 2^/ I . fnn . mnx fnnat\ 
y = ? r— a ^ li sin — sin— r- cos — — . 

Show from this that all the harmonics of the fundamental note of 
the string which correspond to forms of vibration having nodes at 
the point drawn aside by the finger will be wanting in the complex 
note actually sounded. 

Prob. 17.* An iron slab 10 centimeters thick is placed between and 
in contact with two other iron slabs each 10 centimeters thick. The 
temperature of the middle slab is at first 100 degrees Centigrade 
throughout, and of the outside slabs zero throughout. The outer 
faces of the outside slabs are kept at the temperature zero. Re- 
quired the temperature of a point in the middle of the middle slab 
fifteen minutes after the slabs have been placed in contact. 
Given a^ = 0.185 i^ C.G.S. units. Ans. 10^3. 

Prob. 18.* Two iron slabs each 20 centimeters thick, one of which 
is at the temperature zero and the other at 100 degrees Centigrade 
throughout, are placed together face to face, and their outer faces 
are kept at the temperature zero. Find the temperature of a point 
in their common face and of points 10 centimeters from the com- 
:mon face fifteen minutes after the slabs have been put together. 

Ans. 22°.8 ; 15°.! ; i7°.2. 

Art. 12.t Properties of Zonal Harmonics. 

In Art. 4, ^ = Pm{^) was obtained as a particular solution of 
Legendre's Equation [(7), Art. 4] by the device of assuming 
that z could be expressed as a sum or a series of terms of 
-.the form ' a^x"* and then determining the coefficients. We 

♦ See Art. 3. 

f The student should review Art. 4 before beginning this article. 
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PROPERTIES OF ZONAL HARMONICS. 41 

can, however, obtain a particular solution of Legendre's equa- 
tion by an entirely different method. 

The potential function for any point (jr, j, js) due to a unit 
of mass concentrated at a given point (^,,^,, ^i) is 

and this must be a particular solution of Laplace's Equation 
[(3), Art. i], as is easily verified by direct substitution. 
If we transform (i) to spherical coordinates we get 

V= ^ (2) 

Vr^ — 2rr,[cos d cos 0, + sin 6 sin 6/, cos (0—0,)] + r,* 

as a solution of Laplace's Equation in Spherical Coordinates 
[(5). Art. I]. 

If the given point {x^,J/^, z^ is taken on the axis of JT, as it 
must be in order that (2) may be independent of 0, 0^ =0, and 

V=-—= ' - --. (3) 

Vr^ — 2rr, cos ^ + r, 

is a solution of equation (i), Art. 4. 
Equation (3) can be written 

F=i(i-2i:cos^ + ^)-*; (4) 

and if r is less than r, (1—2— cos 6 + —J can be developed 

into a convergent power series. Let 5^«— ^ be this series, 

I t"*' 
p^ being of course a function of B. Then F=— -5^„— ^ is a 

solution of (i), Art. 4. 

Substituting this value of V in the equation, and remem- 
bering that the result must be identically true, we get after a 
slight reduction 
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but, as we have seen, the substitution oi x = cos 6^ reduces this 
to Legendre's equation [(7), Art. 4]. Hence we infer that the 
coefficient of the mth power of z in the development of 
(i — ixz-^-z^)-^ \z a function of x that will satisfy Legendre's 
equation. 

•(i — 2x2 + zy-^ = [i — ^{^^ — z)\ "*» 

and can be developed by the Binomial Theorem ; the coefficient 
of z"^ is easily picked out, and proves to be precisely the func- 
tion of X which in Art. 4 we have represented by PJ^x\ and 
have called a Surface Zonal Harmonic. 
We have, then, 

(I _ 2XZ^!^)-^ = Plx)^Pix) . Z^P^X) . ^ + P,{X) .*•+ . . . (5> 

if the absolute value of z is less than I. 

If ;ir = I, (5) reduces to 
(l-2^ + ^^)-* = P„(l)+P,(.l).i+P,(l).^«+/>3(l).^»+...; 
but (I -2^ + ^')-* = (l -^)-*= 1+^ + ^' + ^' + . . .; 
hence ^«i(0 = i- (6> 

Any Surface Zonal Harmonic may be obtained from the 
two of next lower orders by the aid of the formula 

{n + l)Pn + r{^) - (2« + 1)XPIX) + 71P^_J,X) = O, . (7> 

which is easily obtained, and is convenient when the numerical 
value of X is given. 

Differentiate (5) with respect to z, and we get 

whence 
or by (5) 

(I - 2XZ + Z^){P,{X) + 2Plx) . Z + iPlx) . 2* . . •) 

+■ {z - x){Plx) + Pix) . ^ + PX^) . ^' + . . . 1 = O. (8> 
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Now (8) is identically true, hence the coefficient of each 
power of z must vanish. Picking out the coefficient of ^ and 
writing it lequal to zero, we have formula (7) above. 

By the aid of (7) a table of Zonal Harmonics is easily com- 
puted since we have PJ^x) = i, and PJi^x) = x. Such a table 
for X = cos 6 is given at the end of this chapter. 

Art. 13. Problems in Zonal Harmonics. 

In any problem on Potential if Fis independent of so 
that we can use the form of Laplace's Equation employed in 
Art. 4, and if the value of V on the axis of -A!' is known, and 

•can be expressed as ^a^r^ or as ^ ~^y we can write out 

cour required solution as 

F=2.„r«P.(cos^) or V = ^'^^^ 

iox since PJS) = i each of these forms reduces to the proper 
value on the axis ; and as we have seen in Art. 4 each of them 
satisfies the reduced form of Laplace's Equation. 

As an example, let us suppose a statical charge of M units 
of electricity placed on a conductor in the form of a thin circu- 
lar disk, and let it be required to find the value of the Poten- 
tial Function at any point in the ** field of force** due to the 
•charge. 

The surface density at a point of the plate at a distance r 

irom its centre is 

M 



(T = 



4^;r ^d' — r" 

7t3T 
and all points of the conductor are at potential . See Pierce's 

Newtonian Potential Function (§ 61). 

The value of the potential function at a point in the axis 

ot the plate at the distance x from the plate can be obtained 

without difficulty by a simple integration, and proves to be 

j^ M ,x^ -^ a' 

r = — cos-^-;— — ^. (i) 

2a x' + a' ^ ' 
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44 HARMONIC FUNCTIONS. 

The second member of (i) is easily developed into a power 
series. 

cos - ' — — ; J 

2a ;•:•+«' 

MVn X , x' x" . x' -1 ., 

Hence 



V = ^[5 ^ "Lp^icos e) + ~ Picas 0) 
aL.2 a ^^ ^ ' 3 ^" '^ ^ 



- J J>,(cos^)+...] (4> 



S 
IS our required solution il r <a and 6^ < -, as is 

F = -[- - _ -^/>. (cos ff) + - .-, P,{cos 0) 

-i^>.(cos6i) + ...]ifr>«. (5) 

The series in (4) and (5) are convergent, since they may be 
obtained from the convergent series (2) and (3) by multiplying 
the terms by a set of quantities no one of which exceeds one 
in absolute value. For it will be shown in the next article that 
P^ (cos 0) always lies between i and — i. 

Prob. 19. Find the value of the Potential Function due to the 

attraction of a material circular ring of small cross-section. 

The value on the axis of the ring can be obtained by a sin^ple 

M 
integration, and is . , . , if M is the mass and c the radius of the 
Vc + r 

ring. At any point in space, if r < ^ 

V = y [/'.(cos 0)-j pV.(cos ^) + ^ 7^(008 «)-...], 
and if r > s 
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Art. 14. Additional Forms. 

{a) We have seen in Art. 12 that Pnt{^) is the coefficient of 
;?*• in the development of (i — 2;r^ + ^')-* in a power series. 

(l - 2jr^ + ^) - * = [I - ^(^' + ^ -•') + ^^-l 
= (l — -sr^O" K^ — -srr-^O "*. 

If we develop (i — -&^')"* and {} — -sr^"**)"* by the Bi- 
nomial Theorem their product will give a development for 
(i — 2xz-\-z^)^^. The coefficient of ^"* is easily picked out 
and reduced, and we get 

PJ^cos 0) = 

1.3.5... (J^fn — i) r n \ I . w 

^ ^ ^ ^ —^ 2 COS ni0 + 2 — 7 r cos (m — 2)9' 

2. 4. 6... 2m L i.(2w— i) ^ ^ 

^^ I .2.(2W— l)(2W — 3) ^ ^^ ^ J ^^ 

If w is odd the parenthesis in (i) ends with the term con- 
taining cos 0-, if m is even, with the term containing cos o, but 
in the latter case the term in question will not be multiplied by 
the factor 2, which is common to all the other terms. 

Since all the coefficients in the second member of (i) are 
positive, P„i{cos 0) has its maximum value when ^ = o, and its 
value then has already been shown in Art. 12 to be unity. 
Obviously, then, its minimum value cannot be less than — i. 

{b) If we integrate the value of PJ^x) given in (11), Art. 4, 
m times in succession with respect to jr, the result will be 

found to differ from / , . -(j^ — i)** by terms in- 

{2fn)\ ^ ^ ^ 

volving lower powers of x than the wth. 

Hence P.(.) = -1-, ^(^ - ,)-. (.) 
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{c) Other forms for PJ^x)y which we give without demon- 
stration, are 

p (^) - (zLi): ^ ' (3) 



P«(;.) = 1 /- ./'^ -— .. (f ) 

^i/ [;ir — l/;r' — I . cos 0]**^* 

(4) and (5) can be verified without difficulty by expanding 
and integrating. 

Art. 15. Development in Terms of Zonal Harmonics. 

Whenever, as in Art. 4, we have the value of the Potential 
Function given on the surface of a sphere, and this value de- 
pends only on the distance from the extremity of a diameter, 
it becomes necessary to develop a function of 6 into a series 
of the form 

A^PScos 0) + A,P,(cos 0) + ^,P/cos ^ + . . . ; 

or, what amounts to the same thing, to develop a function of 
X into a" series of the form 

A,P,{x) + A,P,{x) + A^Plx) + .... 

The problem is entirely analogous to that of development 
in sine-series treated at length in Art. 6, and may be solved by 
the same method. 

Assume A^) = A.P^x) + A^P^x) + A^PJ^x) + . . . (i) 

for — I < JT < I. Multiply (i) by PJ^x)dx and integrate from 
— I to I. We get 

1 1 

ff{x)PJ,x)dx J'^lA,fPJx)PJx)dxl f2) 
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We shall show in the next article that 

/ Pm{^)Pn{^)dx = o, unless m^n^ 
-1 

and that /[^«W]V^ = 3-^. 

-1 

Hence A^^""-^^ fAx)P„{x)dx. ' (3) 

-1 

It is important to notice here, as in Art. 6, that the method 
we have used in obtaining A^ amounts essentially to deter- 
mining Af^f so that the equation 

Ax) = A^Plx) + A ,Plx) + AJ>lx) + . . . + A^P„{x) 

shall hold good ior n -{- i equidistant values of x between — i 
and I, and taking its limiting value as n is indefinitely in- 
creased. 

Art. 16. Formulas for Development. 
We have seen in Art. '4 that 2 = PJx) is a solution of 

Legendre's Equation ^ (i — x') —- \-\-m{m -\- i)z = o. (i) 
Hence -^[{i - x') -^- \ + fn{m + i)P^{x) = o, (2) 

Multiply (2) by Pn(^) and (3) by -P^(^), subtract, transpose, 
and integrate. We have 

1 
\tn{m+ I) - «(« + i)-]f PJ,x)P,{x)dx 



=Aw|['--^''^> 
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-1 



-1 
by integration by parts, 

= o. 

1 



Hence f P„{x)P^{x)dx = O, (6) 



-1 

unless m ■=. n. 



If in (4) we integrate from x to i instead of from — i to i, 
we get an important formula. 

(.-.■)[/>.W^-^^-/>.W^] 

PAx)PJ(x)dx^ ^^^- , -X — 7 , X > (7) 

and as a special case, since PJ^^) = i. 

1 (i — x^) 2 
CPSx)dx = -, — ^ , (8) 



unless /« = o. 

1 



To get f[PJix)Ydx is not particularly difficult. By (2), 

-1 
Art. 14, 

By successive integrations by parts, noting that 
T ar^i^* — i)*" contains {x^ — i)* as a factor if /f < w, and 
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that — ^ ,^ ' = {2m)l we get 



/[P^wjv. = (=4,^2?-^ /V - 1)"'^^. (io> 

-1 V V -1 

1 1 

f{x' - i^dx = J'ix - iy\x + i) V4: 
-1 -1 

/« -j- I e/ 

Hehce /[/'„(*)] V^ = ^^^. (i 

-1 ' 

1 

Prob. 2o. Show that / PJ^3c)dx = o if /« is even and is not zero 




m-l 



= (- i)^ -/— r— T . ^'l'^';'"^ , if m is odd. 
' tn\m-\-i) 2.4.6 ... (w — i) 

1 

y» I 
ri'^(^)lVjc = i — , Note that 
"■ ^ ' ■* 2tn -|- I 


\Pfn{po)^ is an even function of x, 

Prob. 22. Show that if f(x) = o from :r = — i to a: = o, and 
f(pc) = I from ^ =• o to :v = i, 

f(x) = i + ^P,(^) - 1-. -P,{x) + i!- . i^ F^x) -... 
2 ' 4 '^ ' 82 12 2.4 '^ ' 

Prob. 23. Show that /?(0) — ^ B„P„{cos S) where 



^„ = ^ULIJ. fp{e)FJcos 0) sinddO. 
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Prob. 24. Show that 

-CSC ^ = y[^ + 5(7)'^. (cos ff) + 9(^-;^) A(cos ff) + .. .]. 
See (i). Art. 14. 

Prob. 25. Show that 



^" = 



1 1 

Note that fx''Pj^x)dx = — i— /'^^ "^K"! '^" ^>y, and use the 

-1 -1 

method of integration by parts freely. 

Prob. 26. Show that if V is the value of the Potential Function 
at any point in a field of force, not imbedded in attracting or repel- 
ling matter; and i( V = f{9) when r =. a, 

V ^ :SaJ^P^{cos e){i r <a 
and V = 2;^^-;;^/'^(cos 6^) if r > a^ 

where A^ = ^^^y/(^)^«(cos 6) sin e^T^. 


Prob. 27. Show that if 

V=c when /•=a; F=^ifr<a, and V= — if r > «. 

r 

Art. 17. Formulas in Zonal Harmonics. 

The following formulas which we give without demonstra- 
tion may be found useful for reference : 

^^=(2«-l)P,..W+(2«-5)/'„_.W+(2«-9)/'„..(4:)+....(l) 

1 

fpixyx = — l-jcp^.x^) - p^^ix)\ (3) 



Digitized by 



Google 



SPHERICAL HARMONICS. 51 

Art. 18. Spherical Harmonics. 

In problems in Potential where the value of Fis given on the- 
surface of a sphere, but is not independent of the angle 0, we 
have to solve Laplace's Equation in the form (5), Art. i, and 
by a treatment analogous to that given in Art. 4 it can be 
proved that 

K = r"* cos nd) sin" ii—-, and k = r*" sin neb sin" 6 — ^ 

where // = cos 6, are particular solutions of (5), Art. i. 

The factors multiplied by r*^ in these values are known as- 
Tesseral Harmonics. They are functions of and 6, and they 
play nearly the same part in unsymmetrical problems that the 
Zonal Harmonics play in those independent of 0. 

F,„(/^, 0) = A^PJ^fA) + "J {A, cos nct> + B, sin «0)sin" o"^^^ 

is known as a Surface Spherical Harmonic of the wth degree^ 
and F=r-F^(/., 0) and F=^F«(yu, 0) 

satisfy Laplace's Equation, (5), Art. I. 

The Tesseral and the Zonal Harmonics are special cases of 
the Spherical Harmonic, as is also a form P,„(cos y) known as 
a Laplace's Coefficient or a Laplacian ; y standing for the angle 
between r and the radius vector r^ of some fixed point. 

For the properties and uses of Spherical Harmonics we 
refer the student to more extended treatises, namely, to 
Ferrer's Spherical Harmonics, to Heine's Kugelfunctionen, or 
to Byerly's Fourier's Series and Spherical Harmonics. 

Art. IQ.*^ Bessel's Funxtions. Properties. 
We have seen in Art. 5 that z = JJ^x) where 



* 



The student should review Art. 5 before reading this article. 
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is a solution of the equation 



(Pz , I dz , , . 

^'+^^+^=°' <^> 

and we have called JJix) a Bessel's Function or Cylindrical 
Harmonic of the zero order. 

,.„=_-^)=£,__fi+__£;____^^+...j(3) 

is called a BesseFs Function of the first order, and 
is a solution of the equation 

which is the result of differentiating (2) with respect to x. 

A table giving values of JJ^x) and /,{x) will be found at 
the end of this chapter. 

If we write JJ^x) for z in equation (2), then multiply 
through by xdx and integrate from zero to x, simplifying the 
resulting equation by integration by parts, we get 

xdJSx) 



dx 



+ fxJlx)dx = o. 



or, since /,(;r) = -^l, 
dx 

X 

fxJlx)dx = xJSx). (5) 



If we write /,(;r) for z in equation (2), then multiply through 
by ^^—j--^ and integrate from zero to x, simplifying by inte- 
gration by parts, we get 

— ' A 



01 



X 

fx{jix)ydx = f [(/„(^))' + {Jix)y\. 



(6) 
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If we replace ;ir by /ir in (2) it becomes 

(See (8). Art. 5). Hence z = JIjjlx) is a solution of (7). 

If we substitute in turn in (j)JJi^Kx) and JJ^^^x) for z, mul- 
tiply the first equation by xJJ^^^x), the second by xjj^^^x), 
subtract the second from the first, simplify by integration by 
parts, and reduce, we get 



JxJ,{^^x)Jlix,x)dx 



Hence if /i^ and /^^ are different roots of JJ^f^a) = o, or of 
/,{^a) = o, or of i^a/,{fAd) — XJS^d) = o, 

a 

J xJXl^Kx)JX^.x)dx == o. (9) 


We give without demonstration the following formulas, 
which are sometimes useful : 

IT 

JJ^x) = - / cos(;ir cos 0)^0. (lo) 



IT 

J^x) = - I sin' cos {x cos 0)^0. (i i) 



They can be confirmed by developing cos {x cos 0), inte- 
grating, and comparing with (i) and (3). 

Art. 20. Applications of Bessel's Functions. 

{d) The problem of Art. 5 is a special case of the following : 
The convex surface and one base of a cylinder of radius a 
and length b are kept at the constant temperature zero, the 
temperature at each point of the other base is a given function 
of the distance of the point from the center of the base ; re- 
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quired the temperature of any point of the cylinder after the 
permanent temperatures have been established. 

Here we have to solve Laplace's Equation in the form 

(see Art, 5), subject to the conditions 

u =zO when ^ = o, 

u =■ o " r = a^ 

u =/{r) " 2 = d. 
Starting with the particular solution of (i), 

u = sinh (/^)/o(/^r), (2) 

and proceeding as in Art. 5, we get, if /i, , yw,, /i,, . . . are roots 
of /oi^a) = o, (3) 

and /(r) = AJ.iM.r) + AJ,{f4,r) + AJlfx/) + . . . , (4) 

(^) If instead of keeping the convex surface of the cylinder 
at temperature zero we surround it by a jacket impervious to 
heat the equation of condition, « = o when r = ^, will be re- 

placed by — = O when r = ^, or if u •=. sinh {jiz)JJ^r) by 

^^ = whenr = «, 

dr 

that is, by — ^/il^a) = o 

or /,{Ma) = o. (6) 

If now in (4) and (5) /^^ , /i,, yu,, . . . are roots of (6), (5) will 
be the solution of our new problem. 

(c) If instead of keeping the convex surface of the cylinder 
at the temperature zero we allow it to cool in air which is at 
the temperature zero, the condition u = when r = a will be 

replaced by ^i^ — [- Au = o when r = a^ A being the coefficient 

of surface conductivity. 
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li u = sinh {m^)/o{m^) this condition becomes 

— M/^{Mr) + A/o{Mr) = o Avhen r = a, 
or M^/iiM^) — ^A/o(/^) = o- (7) 

If now in (4) and (5) /^, , /^, , /^s » • • • ^ire roots of (7), (5) will 
be the solution of our present problem. 

It can be shown that 

/M) = o, (8) 

/M = o, (9) 

and x/,{x) - Uo{^) = O (10) 

have each an infinite number of real positive roots,* The 
earlier roots of these equations can be obtained without serious 
difficulty from the table iov JJ^x) and/,(;r) at the end of this 
chapter. 

Art. 21. Development in Terms of Bessel's Functions. 

We shall now obtain the developments called for in the last 
article. 

Let Ar} = AJJ,lx,r)-^AJSM/) + AJ\{^^^r) + ... (i) 

/^ly M^y Mzy etc., being roots of /Xm^) = o, or of /,(/i^) = o, or 
of M^/iiM^) — Vo(/^^) = o. 

To determine any coefficient Ak multiply (i) by rJJ^^i,r)dr 
and integrate from zero to a. The first member will become 

frf{r)Jln,r)dr. 



Every term of the second member will vanish by (9), Art. 
19, except the term 

a 


frU,{l^>.r)ydr = ±^JxUlx)Jdx- = J ([7.0«*«)]'+|y,(;"*«)]0 



by (6), Art. 19. 



* See Riemann's Partielle Differentialgleichungen, § 97. 
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Hence Aj^ = 



a 

"TTTTT n.^i r /•/ ^\ f^A^Vo(Mkr)dr. (2) 

The development (i) holds good from r = o to r= a (see 
Arts. 6 and 15). 

If Mi9 M^y M99 ^tc, are roots o(/^{/ia) = o, (2) reduces to 

a 

If Mif Mtf M»> etc., are roots o( /^{^a) = o, (2) reduces to 

n 

It /'i* /'.J /^.> etc., are roots of ^J,{^) — 'kj^ixa) = o, 
(2) reduces to 

For the important case where /(r) = i 
frAr)fl^,r)dr=frJlti,r)dr=y^JxJXx)dx = -^/.(/i*«) (6) 





by (5), Art. 19; and (3) reduces to 

2 



(4) reduces to 


A', =0, 


(7) 
(8) 


except for ^ = I 


when pik = 0, and we have 






A, = i; 


(9) 


(5) reduces to 




(10) 



Prob. 28. A cylinder of radius one meter and altitude one meter 
has its upper surface kept at the temperature 100°, and its base and 
convex surface at the temperature 15°, until the stationary temper- 
atures are established. Find the temperature at points on the axis 
25, 50, and 75 centimeters from the base, and also at a point 25 
^centimeters from the base and 50 centimeters from the axis. 

Ans. 29°.6; 47°.6 ; 7i°.2 ; 25°.8 
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Prob. 29. An iron cylinder one meter long and 20 centimeters 
in diameter has its convex surface covered with a so-called non-con- 
ducting cement one centimeter thick. One end and the convex 
surface of the cylinder thus coated are kept at the temperature zero, 
the other end at the temperature of 100 degrees. Given that the con- 
ductivity of iron is 0.185 and of cement 0.00.0162 in C. G. S. units. 

Find to the nearest tenth of a degree the temperature of the mid- 
dle point of the axis, and of the points of the axis 20 centimeters 
from each end after the temperatures have ceased to change. 

Find also the temperature of a point on the surface midway be- 
tween the ends, and of points of the surface 20 centimeters from 
each end. Find the temperatures of the three points of the axis, 
supposing the coating a perfect non-conductor, and again, suppos- 
ing the coating absent. Neglect the curvature of the coating. Ans. 
i5°.4; 4o°.85; 72°.8 ; i5°.3 ; 40^7 ; 72^5 ; o°.o ; o°.o ; i°.3. 

Prob. 30. If the temperature at any point in an infinitely long 
cylinder of radius c is initially a function of the distance of the 
point from the axis, the temperature at any time must satisfy the 

3^ a I'd^u , I 9^\ / . X . . . , , . 

equation ^ = ^ (^ -\ — ^) (see Art. i), since it is clearly in- 
dependent of z and 0. 
Show that 

where, if the surface of the cylinder is kept at the temperature 
zero, /<j , y", , /^, , . . . are roots of J^f^c) = o and Ak is the value 
given in (3) with c written in place of a ; if the surface of the cylin- 
der is adiabatic /i,, /i,, /^j, . . . are roots of /^(mc) = o and A^ is ob- 
tained from (4); and if heat escapes at the surface into air at the tem- 
perature zero A*,, yw,, )W„...are roots of ^cJX^c) — ^/^{^c) = o, 
and Ak is obtained from (5). 

Prob. 31. If the cylinder described in problem 29 is very long 
and is initially at the temperature 100° throughout, and the con- 
vex surface is kept at the temperature 0°, find the temperature of a 
point 5 centimeters from the axis 15 minutes after cooling has begun ; 
iirst when the cylinder is coated, and second, when the coating is 
absent. Ans. 97^.2 ; o°.oi. 

Prob. 32. A circular drumhead of radius a is initially slightly 
distorted into a given form which is a surface of revolution about 
the axis of the drum, and is then allowed to vibrate, and z is the 
ordinate of any point of the membrane at any time. Assuming that 
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2^ must satisfy the equation o7i.= ^ l^ "^~ "" a J» ^^^J^^' ^^ ^'^^ ^^^* 

ditions a: = o when r ■= a^ ^ = o when / = o, and 2: = /(r) when 

<" = o, show that z = AJ^^^) cos /i,^/ + A^jXl^^^) cos yw// + . . . 
where /^j, /^,, /^,, . . . are roots of jXi^) = ^ ^^^^ ^* l^^s the value 
given in (3). 

Prob. 33. Show that if a drumhead be initially distorted as in 
problem 32 it will not in general give a musical note ; that it may be 
initially distorted so as to give a musical note ; that in this case the 
vibration will be a steady vibration ; that the periods of the various 
musical notes that can be given are proportional to the roots of 
J^x) = o, and that the possible nodal lines for such vibrations 
are concentric circles whose radii are proportional to the roots of 

/oW = O. 

Art. 22. Problems in Bessel's Functions. 

If in a problem on the stationary temperatures of a cylinder 
u = o when z = Oy u = o when z = b, and ti =zf[z) when r — a^ 
the problem is easily solved. If in (2), Art. 20, and in the cor- 
responding solution z = cosh {h^)/o{m^) we replace ^ by /Jt, we 
can readily obtain z = sin (jbiz)/Xjiiri) and z = cos (M^)/XMri} 
as particular solutions of (i), Art. 20; and 

and is reaL 

A 2 /»^, . knz _ 
where Aje=^-r f f\z) sm -y- dz f 2> 


by Art. 9. 



r k7rri\ 



I knrix 



Hence « = ^ -^^ sin — 7- ; r- U» 

is the required solution. 
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A table giving the values ol /Jxz) will be found at the end 
of this chapter. 

Prob. 34. A cylinder two feet long and two feet in diameter has 
its bases kept at the temperature zero and its convex surface at 
100 degrees Centigrade until the internabtemperatures have ceased 
to change. Find the temperature of a point on the axis half way 
between the bases, and of a point six inches from the axis, half way 
between the bases. Ans. 72.^1; 8o°.i. 

Art. 23. Bessel's Functions of Higher Order. 

If we are dealing with Laplace's Equation in Cylindrical 
Coordinates and the problem is not symmetrical about an 
axis, functions of the form 

/nW - -^q-^^^ZfTyL . 2\n +iy 2\ 2 \(n + i){n +2) " 'J 
play very much the same part as that played hy JJ{x) in the 
preceding articles. They are known as Bessers Functions of 
the «th order. In problems concerning hollow cylinders much 
more complicated functions enter, known as Bessel's Functions 
of the second kind. 

For a very brief discussion of these functions the reader is 
referred to Byerly's Fourier's Series and Spherical Harmonics ; 
for a much more complete treatment to Gray and Matthews* 
admirable treatise on BesseFs Functions. 

Art. 24. Lamp's Functions. 

Complicated problems in Potential and in allied subjects are 
oisually handled by the aid of various forms of curvilinear co- 
ordinates, and each form has its appropriate Harmonic Func- 
tions, which are usually extremely complicated. For instance. 
Lamp's Functions or Ellipsoidal Harmonics are used when 
solutions of Laplace's Equation in Ellipsoidal coordinates are 
required ; Toroidal Harmonics when solutions of Laplace's 
Equation in Toroidal coordinates are needed. 

For a brief introduction to the theory of these functions 
see Byerly's Fourier's Series and Spherical Harmonics. 
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Table I. Surface Zonal Harmonics. 




B 


P, (cos B) 


P^iCOBB) 


P, (cos B) 


P4 (cos B) 


P5 (cos B) 


P« (cos B) 


Pf (cos B) 


0^ 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1 


.9998 


.9995 


.9991 


.9985 


.9977 


.9967 


.9955 


2 


.9994 


.9982 


.9963 


.9939 


.9909 


.9872 


.9829 


3 


.9986 


.9959 


.9918 


.9863 


.9795 


.9713 


.9617 


4 


.9976 


.9927 


.9854 


.9758 


.9638 


.9495 


.9329- 


5 


.9962 


.9886 


.9773 


.9623 


.9437 


.9216 


.8961 


6 


.9945 


.9836 


.9674 


.9459 


.9194 


.8881 


.8522^ 


7 


.9925 


.9777 


.9557 


.9267 


.8911 


.8476 


.7986. 


8 


.9903 


.9709 


.9423 


.9048 


.8589 


.8053 


.7448. 


9 


.9877 


.9633 


.9273 


.8803 


.8232 


.7571 


.6831 


10 


.9848 


.9548 


.9106 


.8532 


.7840 


.7045 


.6164 


11 


.9816 


.9454 


.8923 


.8238 


.7417 


.6483 


.5461 


12 


.9781 


V .9352 


.8724 


.7920 


.6966 


.5892 


.4732 


13 


.9744 


.9241 


.8511 


.7582 


.6489 


.5273 


.3940 


14 


.9703 


.9122 


.8283 


.7224 


.5990 


.4635 


.3219 


15 


.9659 


.8995 


.8042 


.6847 


.5471 


.3982 


.2454 


16 


.9613 


.8860 


.7787 


.6454 


.4937 


.3322 


.1699' 


17 


.9563 


.8718 


.7519 


.6046 


.4391 


.2660 


.0961 


18 


.9511 


.8568 


.7240 


.5624 • 


.3836 


.2002 


.0289' 


19 


.9455 


.8410 


.6950 


.5192 


.3276 


.1347 


-.0443 


20 


.9397 


.8245 


.6649 


.4750 


.2715 


.0719 


-.1072^ 


21 


.9336 


.8074 


.6338 


.4300 


.2156 


.0107 


-.1662- 


22 


.9272 


.7895 


.6019 


.3845 


.1602 


-.0481 


-.2201 


23 


.9205 


.7710 


.5692 


.3386 


.1057 


-.1038 


-.2681 


24 


.9135 


.7518 


.5357 


.2926 


.0525 


-.1559 


-.3095- 


25 


.9063 


.7321 


.5016 


.2465 


.0009 


-.2053 


-.3463 


26 


.8988 


.7117 


.4670 


.2007 


-.0489 


-.2478 


-.3717 


27 


.8910 


.6908 


.4319 


.1553 


-.0964 


-.2869 


-.3921 


28 


.8829 


.6694 


.3964 


.1105 


-.1415 


-.3211 


-.4052 


29 


.8746 


.6474 


.3607 


.0665 


-.1839 


-.3503 


-.4114 


30 


.8660 


.6250 


.3248 


.0234 


-.2233 


-.3740 


-.4101 


31 


.8572 


.6021 


.2887 


-.0185 


-.2595 


-.3924 


-.4022- 


32 


.8480 


.5788 


.2527 


-.0591 


-.2923 


-.4052 


-.3876- 


33 


.8387 


.5551 


.2167 


-.0982 


-.3216 


-.4126 


-.3670- 


34 


.8290 


.5310 


.1809 


-.1357 


-.3473 


-.4148 


-.3409- 


35 


.8192 


.5065 


.1454 


-.1714 


-.3691 


-.4115 


-.3096 


36 


.8090 


.4818 


.1102 


-.2052 


-.3871 


-.4031 


-.2788 


37 


.7986 


.4567 


.0755 


-.2370 


-.4011 


-.3898 


-.2343 


38 


.7880 


.4314 


.0413 


-.2666 


-.4112 


-.3719 


-.1918 


39 


.7771 


.4059 


.0077 


-.2940 


-.4174 


-.3497 


-.146^ 


40 


.7660 


.3802 


-.0252 


-.3190 


-.4197 


-.3234 


-.1003 


41 


.7547 


.3544 


-.0574 


-.3416 


-.4181 


-.2938 


-.0534 


42 


.7431 


.3284 


-.0887 


-.3616 


-.4128 


-.2611 


-.0065 


43 


.7314 


.3023 


-.1191 


-.3791 


-.4038 


-.2255 


.0398 


44 


.7193 


.2762 


-.1485 


-.3940 


-.3914 


-.1878 


.0846 


45" 


, .7071 


.2500 


-.1768 


-.4062 


-.3757 


-.1485 


.1270 
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Table I. Surface Zonal Harmonics. 




B 


P, (cos B) 


Pa cos B) 


P, (cos B) 


P4(C08tf) 


P5 (cos B) 


P.(C08«) 


PT(COg0) 


45^ 


.7071 


.2500 


-.1768 


-.4062 


-.3757 


-.1485 


.1270 


46 


.6947 


.2238 


-.2040 


-.4158 


-.3568 


-.1079 


.1666 


47 


.6820 


.1977 


-.2300 


-.4252 


-.3350 


-.0645 


.2054 


48 


.6691 


.1716 


-.2547 


-.4270 


-.3105 


-.0251 


.2349 


49 


.6561 


.1456 


-.2781 


-.4286 


-.2836 


.0161 


.2627 


50 


.6428 


.1198 


-.8002 


-.4275 


-.2545 


.0563 


.2854 


51 


.6293 


.0941 


-.8209 


-.4239 


-.2235 


.0954 


.3031 


52 


.6157 


.0686 


-.3401 


-.4178 


-.1910 


.1326 


.3153 


53 


.6018 


.0433 


-.3578 


-.4093 


-.1571 


.1677 


.8221 


54 


.6878 


.0182 


-.3740 


-.3984 


-.1223 


.2002 


.3284 


55 


.5736 


-.0065 


-.3886 


-.3852 


-.0868 


.2297 


.3191 


56 


.5592 


-.0310 


-.4016 


-.3698 


-.0510 


.2559 


.3095- 


57 


.5446 


-.0551 


-.4131 


-.3524 


-.0150 


.2787 


.2949 


58 


.5299 


-.07>8 


-.4229 


-.3331 


.0206 


.2976 


.2752 


59 


.5150 


-.1021 


-.4310 


-.3119 


.0557 


.3125 


.2511 


60 


.5000 


-.1250 


-.4375 


-.2891 


.0898 


.3232 


.2231 


61 


.4848 


-.1474 


-.4423 


-.2647 


.1229 


.3298 


.1916 


62 


.4695 


-.1694 


-.4455 


-.2390 


.1545 


.3321 


.1571 


63 


.4540 


-.1908 


-.4471 


-.2121 


.1844 


.3302 


.1203 


64 


.4384 


-.2117 


-.4470 


-.1841 


.2123 


.3240 


.0818 


65 


.4226 


-.2821 


-.4452 


-.1552 


.2381 


.3138 


.0422 


66 


.4067 


-.2518 


-.4419 


-.1256 


.2615 


.2996 


.0021 


67 


.3907 


-.2710 


-.4370 


-.0955 


.2824 


.2819 


-.0375 


68 


.3746 


-.2896 


-.4305 


-.0650 


.3005 


.2605 


-.0763 


69 


.3584 


-.3074 


-.4225 


-.0344 


.3158 


.2361 


-.1135 


70 


.3420 


-.3245 


-.4130 


-.0038 


.3281 


.2089 


-.1485 


71 


.3256 


-.3410 


-.4021 


.0267 


.3373 


.1786 


-.1811 


72 


.3090 


-.3568 


-.3898 


.0568 


.3434 


.1472 


-.2099 


73 


.2924 


-.8718 


-.3761 


.0864 


.3463 


.1144 


-.2347 


74 


.2756 


-.3860 


-.3611 


.1153 


.3461 


.0795 


-.2559 


75 


.2588 


-.3995 


-.3449 


.1434 


.3427 


.0481 


-.2730 


76 


.2419 


-.4112 


-.3275 


.1705 


.3362 


.0076 


-.2848 


77 


.2250 


-.4241 


-.3090 


.1964 


.3267 


-.0284 


-.2919 


78 


.2079 


-.4352 


-.2894 


.2211 


.8143 


-.0644 


-.2943 


79 


.1908 


-.4454 


-.2688 


.2448 


.2990 


-.0989 


-.2918 


80 


.1736 


-.4548 


-.2474 


.2659 


.2810 


-.1321 


-.2835 


81 


.1564 


-.4633 


-.2251 


.2859 


.2606 


-.1635 


-.2709 


82 


.1392 


-.4709 


-.2020 


.3040 


.2378 


-.1926 


-.2536 


83 


.1219 


-.4777 


-.1783 


.3203 


.2129 


-.2193 


-.2321 


84 


.1045 


-.4836 


-.1589 


.8345 


.1861 


-.2431 


-.2067 


85 


.0872 


-.4886 


-.1291 


.3468 


.1577 


-.2638 


-.1779 


86 


.0698 


-.4927 


-.1038 


.3569 


.1278 


-.2811 


-.1460 


87 


.0523 


-.4959 


-.0781 


.3648 


.0969 


-.2947 


-.1117 


88 


.0849 


-.4982 


-.0522 


.3704 


.0651 


-.3045 


-.0735 


89 


.0175 


-.4995 


-.0262 


.3739 


.0327 


-.3105 


-.0381 


90° 


.0000 


-.5000 


.0000 


.3750 


.0000 


-.3125 


.0000 
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Table II. Bbssel's Functions. 



X 


Joix) 


Jiix) 


X 


Joix) 


^ Jiix) 


X 


^o(^) 


Jiix) 


0.0 


1.0000 


0.0000 


5.0 


-.1776 


-.3276 ; 


10.0 


-.2459 


.0435 


0.1 


.9975 


.0499 


5.1 


-.1443 


-.3371 ' 


10.1 


-.2490 


.0184 


0.2 


.9900 


.0995 i 


5.2 


-.1103 


-.3432 


10.2 


-.2496 


.0066 


0.8 


.9776 


.1483 


5.3 


-.0758 


-.3460 


10.3 


-.2477 


-.0313 


0.4 


.9604 


.1960 


5.4 


-.0412 


-.3453 


10.4 


-.2434 


-.0555 


0.5 


.9385 


.2423 


5.5 


-.0068 


-.3414 


10.5 


-.2366 


-.0789 


0.6 


.9120 


.2867 


5.6 


.0270 


-.3343 


10.6 


-.2276 


-.1012 


0.7 


.8812 


.3290 


5.7 


.0599 


-.8241 


10.7 


-.2164 


-.1224 


0.8 


.8463 


.3688 


5.8 


.0917 


-.3110 


10.8 


-.2032 


-.1422 


0.9 


.8075 


.4060 


5.9 


.1220 


-.2951 


10.9 


-.1881 


-.1604 


1.0 


.7652 


.4401 


6.0 


.1506 


-.2767 


11.0 


-.1712 


-.1768 


1.1 


.7196 


.4709 1 


6.1 


.1773 


-.2559 


11.1 


-.1528 


-.1913 


1.2 


.6711 


.4983 


6.2 


.2017 


-.2329 


11.2 


-.1330 


-.2039 


1.8 


.6201 


.5220 


6.3 


.2238 


-.2081 


11.3 


-.1121 


-.2143 


1.4 


.5669 


.5419 


6.4 


.2433 


-.1816 


11.4 


-.0902 


-.2225 


1.5 


.5118 


.5579 ; 


6.5 


.2601 


-.1538 


11.5 


-.0677 


-.2284 


1.6 


.4554 


.5699 1 


6.6 


.2740 


-.1250 


11.6 


-.0440 


-.2320 


1.7 


.3980 


.5778 


6.7 


.2851 


-.0953 


11.7 


-.0213 


-.2333 


1.8 


.3400 


.5815 


6.8 


.2931 


-.0652 


11.8 


.0020 


-.2323 


1.9 


.2818 


.5812 


6.9 


.2981 


-.0349 


11.9 


.0250 


-.2290 


2.0 


.2239 


.5767 


7.0 


.3001 


-.0047 


12.0 


.0477 


-.S234 


2.1 


.1666 


.5688 1 


7.1 


.2991 


.0252 


1 12.1 


.0697 


-.2157 


2.2 


.1104 


.5560 


7.2 


.2951 


.0543 


1 12.2 


.0908 


-.2(60 


2.3 


.0555 


.6399 


7.3 


.2882 


.0826 


12.3 


.1108 


-.1943 


2.4 


.0025 


.5202 


7.4 


.2786 


.1096 


12.4 


.1296 


-.1807 


2.5 


-.0484 


.4971 


7.5 


.2663 


.1352 


12.5 


.1469 


-.1655 


2.6 


-.0968 


.4708 


7.6 


.2516 


.1592 


12.6 


.1626 


-.1487 


2.7 


-.1424 


.4416 


7.7 


.2346 


.1813 


12.7 


.1766 


-.1307 


2.8 


-.1850 


.4097 


7.8 


.2154 


.2014 


12.8 


.1887 


-.1114 


2.9 


-.2243 


.3754 


7.9 


.1944 


.2192 


12.9 


.1988 


-.0912 


3.0 


-.2601 


.3391 


8.0 


.1717 


.2346 


13.0 


.2069 


-.0703 


3.1 


-.2921 


.3009 


8.1 


.1475 


.2476 


! 13 1 


.2129 


-.0489 


3.2 


-.3202 


.2613 


82 


.1222 


.2580 


13.2 


.2167 


-.0271 


3.3 


-.3443 


.2207 


8.3 


.0960 


.2657 


13.3 


.2183 


-.0052 


3.4 


-.3643 


.1792 


8.4 


.0692 


.2708 


13.4 


.2177 


.0166 


3.5 


-.3801 


.1874 


8.5 


.0419 


.2731 


13.5 


.2150 


.0380 


8.6 


-.3918 


.0955 


8.6 


.0146 


.2728 


13.6 


.2101 


.0590 


3.7 


-.3992 


.0538 


8.7 


-.0125 


.2697 


13.7 


.2032 


.0791 


3.8 


-.4026 


.0128 


8.8 


-.0392 


.2641 


13.8 


.1943 


.0984 


3.9 


-.4018 


-.0272 


8.9 


-.0653 


.2559 


13.9 


.1836 


.1166 


4.0 


-.3972 


-.0660 


9.0 


-.0903 


.2453 


14.0 


.1711 


.1384 


4.1 


-.3887 


-.1033 


9.1 


-.1142 


.2324 


14.1 


.1570 


.1488 


4.2 


-.3766 


-.1386 


9.2 


-.1367 


.2174 


14.2 


.1414 


.1626 


4.3 


-.3610 


-.1719 


9.3 


-.1577 


.2004 


14.3 


.1245 


.1747 


4.4 


-.3423 


-.2028 


9.4 


-.1768 


.1816 


14.4 


.1065 


.1850 


4.5 


-.3205 


-.2311 


9.5 


-.1939 


.1613 


14.5 


.0875 


.1934 


4.6 


-.2961 


-.2566 


9.6 


- .2090 


.1395 


14.6 


.0679 


.1999 


4.7 


-.2693 


-.2791 


9.7 


-.2218 


.1166 


14.7 


.0476 


.2043 


4.8 


-.2404 


-.2985 


9.8 


-.2823 


.0928 


14.8 


.0271 


.2066 


4.9 


-.2097 


-.3147 


9.9 


-.2403 


.0684 


14.9 


.0064 


.20fi9 


5.0 


-.1776 


-.8276 


10.0 


-.2459 


.0435 


15.0 


-.0142 


.2051 
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n 


x^ for JoixJ = 


x^ for JiixJ = 


n 


x^ for JoixJ = 


x,forJi(a;^ = fl 


1 
2 
3 
4 
5 


2.4048 

6.5201 

8.6537 

11.7915 

14.9309 


3 8817 

7.0156 

10.1785 

18.3287 

16.4706 


6 

7 

8 

9 

10 


18.0711 
21.2116 
24.3625 
27.4936 
30.6346 
• 


19.6169 
22.7601 
26 9087 
29.0468 
32.1897 



Table IV.— Values op Jo{xi), 



X 


Joixi) 


x< 


Joixi) 


X 


Joixi) 


0.0 


1.0000 


2.0 


2 2796 


4.0 


11.8019 


0.1 


1.0025 


2.1 


2.4468 


4.1 


12.3286 


0.2 


1.0100 


2.2 


2.6291 


4.2 


18.4425 


0.3 


1.0226 


2.3 


2.8296 


4.8 


14.6680 


04 


1.0404 


2.4 


8.0493 


4.4 


16.0104 


0.5 


1.0635 


2.5 


3.2898 


4.6 


17.4812 


0.6 


1.0920 


2.6 


3.6533 


4.6 


19.0926 


7 


1.1263 


2.7 


3.8417 


4.7 


20.8685 


0.8 


1.1665 


2.8 


4 1573 


4.8 


22.7937 


0.9 


1.2130 


2.9 


4.5027 


4.9 


24.9148 


1.0 


1.2661 


3.0 


4.8808 


6.0 


27.2899 


1.1 


1.3262 


3.1 


6.2945 


6.1 


29.7889 


1.2 


1.3937 


3.2 


6.7472 


5.2 


32.5886 


1.3 


1.4968 


3.3 


6.2426 


6.8 


35.6481 


14 


1.5534 


8.4 


6.7848 


5.4 


89.0088 


1.6 


1.6467 


8.5 


7.3782 


6.6 


42.6946 


1.6 


1.7600 


3.6 


8.0277 


5.6 


46.7876 


1.7 


1.8640 


8.7 


8.7886 


6.7 


51.1725 


1.8 


1.9896 


3.8 


9.6169 


6.8 


66.0881 


1.9 


2.1277 


3.9 


10.8690 


6.9 


61.3766 



Digitized by 



Google 



Digitized by 



Google 



INDEX. 



Bemouilli, Daniel, 7. 
BessePs Functions: 
applications to physical problems, 

53-55. 
development in terms of, 55-56. 
first used, 7. 

introductory problem, 21. 
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Fourier, 7. 
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63- 
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tesseral, 51. 

toroidal, 59. 
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Legendre, 7. 
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Legendre's equation, 17, 40, 41, 47. 

Musical strings, 7. 
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problems, 39, 40. 

Perry, John, 8. 
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problems, 44, 51. 

Sine series, 26. 
determination of the coefficients, 26- 

28. 
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Stationary temperatures: 
problems, 21, 25, 56, 57, 59. 
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Tables, 60-63. 
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